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This revision of the Floating Point Interpreter User's Man- 
ual, 093-000019-04, supersedes 093-000019-03 and con- 
stitutes a minor revision. A list of changes is given at the 
back of the manual, following the index. 



Introduction 

Most small, general purpose computers do not have hardware for the manipulation 
of floating point numbers. They have, therefore, implemented software packages 
to provide floating point functions. Two approaches can be taken to solve this soft- 
ware problem. 

The first approach is to provide a number of subroutines, each of which performs 
a specific function, e.g. floating addition. This requires that the user pass inputs 
in a standard manner to the subroutine. Since many small computers have only 
two accumulators, only a single operand (represented in two words) can be passed 
in the accumulators. The address of the second operand must be provided with 
the subroutine call. A typical calling sequence requires two load accumulator 
instructions, a subroutine call, and an address word. One operand is always 
destroyed and replaced by the result (usually the accumulator operand). This 
means that an intermediate result cannot be tested without destroying an operand. 
A lengthy program requiring these subroutines cannot make efficient use of 
storage and the manipulation of operands becomes a tedious, cumbersome job. 

The second approach is to provide an interpreter in which all the floating point 
functions are imbedded and which can simulate floating point registers in storacrp 
The interpreter is given control by a subroutine call. "6nce in control, it accesses 
succeeding memory words and "interprets" a sixteen bit word in a manner similar 
to the hardware. These instructions are not executed in the hardware sense but 
are instead interpreted to provide extended machine features. The interpreter 
approach was adopted for the Nova- line floating point package. 

The interpreter provides four floating point accumulators which can be addressed 
and manipulated in a manner similar to the actual hardware accumulators. 
Floating point instructions are syntactically similar to machine instructions* and 
are assembled in a similar way. For example, the instruction: 

ADD* 0,2,SNR 

which adds AC0 to AC2, skips the next instruction if the result is non-zero, and 
does not change AC0 or AC2, has a similar floating point version. The instruction 

FADD# 0.2.FSNR 

adds floating accumulator (designated as FAC0) to FAC2, skips the next instruction 
if the result is non-zero, and does not destroy FAC0 or FAC2. 

The interpreter can be implemented in Read Only Storage (ROS), and is completely 
reentrant. These properties, in addition to the floating point instruction set, are 
described in the body of this manual. 



The reader should be thoroughly familar with machine instructions and the 
Nova assembler. Instructions are described in "How to Use the Nova 
Computers"; the absolute and relocatable assemblers are described in 
documents 093-000017 and 093-000040 respectively. 
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STRUCTURE OF THE INTERPRETER 
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The interpreter is non-destructive, i.e., no instruction 
is modified in any manner during execution. Further, it 
does not store temporary information within itself but 
uses a writable area which must be provided by the user. 

Within the interpreter, subroutine linkage is via a push- 
down list, maintained in the writable area provided by 
the user. This property, coupled with the properties 
already mentioned, makes the interpreter reentrant. 
This means that if a second user routine requires the 
interpreter, it may interrupt the current routine, per- 
form its function using the interpreter, and return to the 
first routine without affecting the state of the interpreter. 

The absolute interpreters (basic binary and extended 
binary) require eight page locations which the user 
must not destroy. These locations and their uses are: 

004 contains the starting address of the interpreter. 

005 contains the starting address of the initializa - 
tion routine for the interpreter. 

006 one temporary word used by the interpreter and 
saved by any routine reentering the package. 

007 a word set up by the user to contain the address 
of a writable area for the interpreter. 

040 a word containing the address of a user written 
'get character" routine. 

041 a word containing the address of a u.ier written 
"put character" routine. 

042-043 words containing the linkage addresses to deter- 
mine the location of extended code (if present). 



Floating point numbers are internally stored in two 
consecutive 16-bit words. The form is: 



1 



K 



1 8 



31 



M 



S is the sign of the mantissa, M, in bits 8-31. The 
mantissa is considered to be a normalized six digit 
hexadecimal fraction, and the range of the magnitude 
of the mantissa is: 

16**-1 «M«(l-16**-6). 

The characteristic, C, is the intege r exponent of 16 
in excess 64i0 code. The total range of magnitudes 



16**-l*16**-63^ F $ (1-16**-6)*16**63 

or approximately 

2.4*10**-78 « F < 7.2*10**75 

Any operand having a zero mantissa is represented in 
true zero form, i.e. , bits 0-31 are 0. Negative 
numbers are identical to their positive counterparts 
except S = 1 instead of 0. 

The maximum error of a normalized mantissa is 
less than 16**-6. 



In the relocatable version, the same functions are pro- 
vided by six 7REL locations and two NREL locations. 
Those of interest to the user are ITMP, which corres- 
ponds to location 6 in the absolute versions and is de- 
clared as an entry by the interpreter, and WSA, GETC 
and PUTC, which correspond to locations 7, 40, and 41 
respectively and which must be page zero relocatable. 
WSA, GETC and PUTC are defined as entries by the user 
program. 



See also "How to Use the NOVA Computers", Appendix C. 
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GENERAL 

The interpreter provides four floating point accumula- 
tors. They are numbered 0, 1, 2, and 3 like the 
hardware accumulators. The designation FACn will 
be used for floating accumulator n. Arithmetic is per- 
formed accumulator-to-accumulator as with fixed 
point instructions. Operands can be accessed and 
stored using memory reference instructions. Instruc- 
tions which reference floating point operands in memory 
should provide an address which points to the first word 
of the two word operand. If indexing is specified, the 
hardware index register is used. For example, 
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1,4,2 



loads FAC1 from two consecutive words in memory 
whose fix si word address is 4 -t- C(AC2). * Certain 
instructions that manipulate hardware accumulators 
AC2 and AC3 will be described in the section 'Floating 

vided for manipulating AC0 and AC1 with the floating 
point instruction set. 

BASIC AND EXTENDED INTERPRETERS 

This chapter describes the features of the Basic Float- 
ing Point Interpreter. An extended version of the in- 
terpreter, in both binary** and relocatable format, is 
available. The extended version has the same features 
described here; additional features of the Extended and 
Extended Relocatable Interpreters are described in 
Chapters 3 and 4 respectively. 

THE WRITABLE AREA 

The Basic Floating Point Interpreter requires 64 (decimal) 
words of contiguous writable memory. If the Extended 
Interpreter is used, 100 (decimal) words must be pro- 
vided. The Extended Interpreter and its added features 
are described in detail in Chapter 3. Other than the 
difference in the length, the writable area of both 
versions is the same. 

The address of the first word of the writable area pro- 
vided" by the user must be stored in location 007 of 
page before any interpreter commands are executed. 
If a second routine (or third, etc. ) may reenter the 
interpreter, the routine must provide an address 
pointing to a different memory area before reentrance is 
made. 



Basic Interpreter Use and Command Structure 

THE WRITABLE AREA (cont'd) 

A number of flags are stored in this writable area and 
may be examined by the user. To access the flags, 
the user must exit from the interpretive mode (des- 
cribed on page 2-2). if an index register is loaded 
» iuu] "A-tiUuu vpi> uii; uoci niciy access Liie nag 
words by an instruction of the form: 

LDA 0,n,i_ 

where: _i is the index register (2 or 3) . 

n is a constant displacement as described 
below. 

n = 0: The first word contains the overflow/under- 
flow flags. If the result of an operation is 
less than 16**-64, bit 15 of this word will be 
set to indicate underflow. If the result of an 
operation is greater than (l-16**-6)*16**+63, 
bit 14 of the word will be set to indicate over- 
flow. Other conditions that may set overflow/ 
underflow will be discussed where appropriate. 
This word is initially cleared. Once a bit is set 
it will remain set. It is the user's responsibil- 
ity to reset these bits. 

n = 1: After an input conversion (see page 2 -6\ this 
word will be zero if n<3 conversion was per- 
formed (because of an input error). Other- 
wise the word will be non-zero. 

n = 2: After an input conversion, this word will con- 
tain the 7-bit ASCII character, right justified, 
that served as the break character in the input 
stream. 

n = 3: After an input conversion, this word will be 
if no decimal point (". ") was encountered in 
the input stream. If a decimal point was seen, 
the word will be 1. 

The diagram below summarizes these flags: 

8 9 1415 word 



007 



COK 



~E3p=£± 



BRKC 



used only by 
interpreter 




OFLO 
UFLO 

FPTF 



C(x) means 'contents of x'. 

The Basic Interpreter (binary) requires locations 56008 " 76 008 » a nd the Extended Interpreter 
(binary) requires locations 41008 



7600b 
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Chapter 2- Basic Interpreter Use anil Command Structure 

THE WRITABLE AREA (cont'd) 

A convenient means of accessing the flags is to define 
four symbolic equivalences such as: 

FLGS = ; OVERFLOW/UNDERFLOW FLAGS 

COK = 1 ; CONVERSION OK 

BRKC=2 ; BREAK CHARACTER 

FPTF= 3 ; FLOATING POINT FLAG 

Now if AC2 contains the address of the writable area, 
any flag may be accessed by statements like the 
following: 

LDA 0, FLGS, 2 

LDA l,COK,2 

LDA 1.BRKC.2 

LDA 1.FPTF.2 

INITIALIZATION 

The interpreter, or more correctly the writable area, 
must be initialized before floating point instructions 
are executed. Initialization should be given once for 
every writable area to be used. The command is 

FINI 



which generates the instruction 

JSR @5 

(Note that the initialization routine address is in loca- 
tion 005. ) Location 007 must point to the writable 
area. This command destroys AC3, but preserves 
all other accumulators and Carry. The initialization 
routine clears the overflow /underflow flags and sets up 
linkage for the push-down list. 

ENTERING-EXITING THE INTERPRETIVE MODE 



To use the Floating Point Interpreter, it is necessary 
to distinguish between the processing of normal instruc- 
tions (described in "How to Use the Nova Computers") 
and the processing of floating point instructions. Since 
floating point instructions are "interpreted" (not 
executed per se), the interpreter must be given control 
before floating point code will be executed properly. 
Whenever the interpreter is in control, this will be 
called the "interpretive mode. " Otherwise, the machine 
will be referred to as in "normal mode. " To enter the 
interpretive mode, the command is 

FETR 



ENTERING-EXITING THE INTERPRETIVE MODI': (cmii'd) 

which generates the instruction: 

JSR @4 

As noted earlier, location 004 will contain the starting 
address of the interpreter. Once in the interpretive 
mode, only the floating point instruction set can be 
used. (Normal instructions will be decoded and exe- 
cuted as floating point instructions. ) To return to 
normal processing, the command 

FEXT 

must be given. 

FETR destroys AC3, but the contents of all other 
accumulators and Carry will be saved. AC2 and AC3 can 
be used for indexing in the interpretive mode. Certain 
instructions also enable modification of the contents of 
AC2 and AC3. FEXT will restore AC0, AC!, and Carry to 
their state before entering the interpreter. AC2 and AC3 
will reflect any changes caused by floating point 
instructions that modified their contents. 

The normal sequence of code using the interpreter would 
be of the form : 

FINI INITIALIZE INTERPRETER 
LDA 0.CNST 

STA 0,TEMP ;NORMAL INSTRUCTIONS 

FETR ;ENTER INTERPRETER 

; FLOATING POINT INSTRUCTIONS 



FEXT 
LDA 1, C5 



;EXIT INTERPRETER 



FLOATING POINT INSTRUCTION SET 

The set of instructions for the Basic Floating Point 
Interpreter obey the syntactic rules described in the 
Assembler Manual, 093-000017. For example, all 
floating ALC instructions require a source accumula- 
tor and a destination accumulator. All floating point 
instructions begin with an "F" to distinguish them 
from normal instructions. Appendix A summarizes 
the instructions and Appendix B gives their octal en - 
coding. 
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Memory Reference Instructions 

All memory reference instructions require an address, 
which is represented in the following discussion by 
" adr ". The effective address is calculated exactly as in 
the standard Nova instruction set, except that the auto- 
incrementing, auto -decrementing properties of locations 
20-37 have not been implemented. (For a discussion of 
effective address calculation, see "How to Use the Nova 
Computers", Section 2-1. ) 

The instruction 

FLDA n, adr 

causes FACn to be loaded with the two word operand at 
adr, adr+ 1. 

The instruction 

FSTA n, adr 
causes FACn to be stored in memory at adr , adr+1. 
The instruction 

FJMP adr 

causes control to be transferred to the floating point 
instruction at adr. 

The instruction 

FJSR adr 

causes control to be transferred to the instruction 

at prlr anrl AT^ tr» Ue> oot *-<-» t-U^ i7nln A «* «-i — *. 

— — — . .. — — w. ^^.^ t ,j ^,^ ut-t tw \.xx\, vaiut ui tile CUX±CilL 

location counter + 1. For example, a floating point 
subroutine can be executed by the following: 



FLDA 0, LOC 
FJSR SUBR 
FSTA 0, RSLT 



;MAIN PROGRAM 

;LOAD FAC0 

JUMP TO SUBROUTINE 

;STORE RESULT 



Chapter 2 - Basic Interpreter Use and Command Structure 

Memory Reference Instructions (cont'd) 

converted to a fixed point, double precision integer 
(truncated) at adr, adr+1 . If the conversion results in 
an integer whose absolute value is greater than 
2**24-1, the overflow flag will be set and 2**24-1 will 
be returned as the magnitude. Note that while floating 
point numbers are represented in signed- magnitude 
format, fixed point values will always be represented 
in two's complement notation. 

The instruction 

FFLOadr 

causes the fixed point, double precision integer at 
adr , adr+ 1 to be converted to a floating point number 
at adr, adr+ 1. Negative integers must be represented 
in two's r.omnlement format. 

The instruction 

FISZ adr 

causes the contents of adr to be incremented by one and 
the next floating point instruction in sequence to be 
skipped if the result is zero. 

The instruction 

FDSZ adr 

causes the contents of adr to be decremented by one 
and the next floating point instruction in sequence to 
be skipped if the result is zero. FISZ and FDSZ should be 
used with fixed point, single precision integers - not 
with floating point numbers. 

The instruction 

FST3 adr 

causes AC3 to be stored at adr . 

The instruction 



SUBR: 



;SUBROUTINE 



FJMP 0, 3 ;RETURN 

The instruction 

FFIXadr 
causes the floating point number at adr, adr+1 to be 



FLD3 adr 

causes AC3 to be loaded from the contents of adr. FST3 
and FLD3 operate on real accumulator 3. FST3 is 
useful for saving the return address inside a floating point 
subroutine. Note that the return address should always 
be saved if a floating point subroutine exits and then 
enters the interpretive mode, since FETR destroys AC3. 
These two instructions also provide a means for initial- 
izing loop counts without leaving the interpretive mode. 
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Chapter 2 - Basic Interpreter Use and Command Structure 

Memory Reference Instructions (cont' d) Arithmetic Instructions (ALC) (cont" cl) 

For example, the sequence: 

FEXT 

LDA 0.CNT 

STA 0.TEMP 

FETR 



LP1: 



FEXT 

DSZ TEMP 
JMP LP1 
FETR 



can be replaced by: 



LP1: 



FLD3 CNT 
FST3 TEMP 



FDSZ TEMP 
FJMP LP1 

Arithmetic Instructions (ALC) 

Options 

The floating point ALC instructions are similar to 
normal ALC instructions. Two floating accumulators 
must be specified. The first is the source accumulator, 
the second the destination accumulator. 

Seven skip conditions are defined (in addition to the 
default "no skip"). These conditions are listed in 
Table 1. The conditions FSZR and FSNR should be 
used with caution. Since floating point arithmetic is 
inner', nth approximate, the probability of obtaining 
true zero is very low. The normal procedure is to 
test a result (or difference) against a small quantity, 
£ . For example, if we wish to test for the con- 
vergence of an iterative procedure we might use the 
following: 



EPSLN: 1. 0E-6 



;EPSILON IS 10** -5 



Skip Mnemonic 


Effect 


FSLT 


skip if result < 


FSLE 


skip if result < 


FSGT 


skip if result > 


FSGE 


skip if result >, 


FSNR 


skip if result f 


FSZR 


skip if result = 


FSKP 


unconditional skip 



Table 1 - Skip Mnemonics 

All ALC instructions permit the load/no load option. 
As with normal instructions, if a floating point instruc- 
tion mnemonic is suffixed with "#", the results of the 
operation will not replace the contents of the destination 
register. 

A further option is available with one class of ALC 
instructions. This option will prevent post-normaliza- 
tion of the result. The instructions described in the 
section following permit this option. It is called for 
by suffixing "U" (for unnormalized) to the instruction 
mnemonic. For example, FMOV 0, 1 moves FAC0 , 
to FAC1 and normalizes the result, while FMOVU 0, 1 
moves FAC0 to FAC1 without normalization. 

ALC Instructions with Post-Normalize Option 

The instruction 

FMOV n,m 
moves FACn to FACm. 
The instruction 

FPOS n, m 
moves the absolute value of FACn to FACm. 



FLDA 


0, ORSLT 


GET OLD RESULT 


FLDA 


1, NRSLT 


GET NEW RESULT 


FSUB 


1, 


OLD- NEW 


FLDA 


1, EPSLN 


EPSILON 


FPOS 


0,0 


ABS (OLD-NEW) 


FSUB 


1,0, FSLE 


SKIP IF < EPSILON 



The instruction 

FMNS n, m 

moves the negative of the absolute value of FACn to 
FACm. 
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ALC Instructions with Post-Normalize Option- (cont'd) 
The instruction 

FNEG n, m 

"*»-» • ^-"-J U1V llCgULtl ¥ V- * U. t l_l\- VJ-L 1 J1\J 11 |,K^ 1 JIVjill, 

The instruction 

FRND n,m 

rounds the value of FACn and moves it to FACm. By 
"round" we mean the following. 



r - iu • ■ -o ' jio ■ o r-t-i/2j 
The instruction 



,<"> 
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ALC Instructions that Always Post-Normalize 

This class of ALC instructions always post-normalizes 
the result. They assume pre-normalized operands. 
Overflow is checked and indicated by setting the over- 
flow flag and returning a magnitude of (1-16**-6)*16**63 
as the magnitude of the result. Underflow is checked 
and indicated by setting the underflow flag and returning 
true as the result. 

The instruction 

FMPY n,m 

multiplies FACn by FACm and moves the result to 
FACm. 

The instruction 



r/ijjjj n. m 



adds FACn to FACm and moves the result to FACm. 
L unuerflow occurs, the underflow flag is set and true 
is returned as the result. If overflow occurs, the 
overflow flag is set and a magnitude of (1-16**-6)*16*H63 
is returned as the result. The operands are assumed 
to be pre-normalized. 

The instruction 



FSUB n,m 

subtracts FACn from FACm and moves the result to 
FACm. The overflow conditions are handled as 
with FADD. Prenormalized operands are assumed. 

Table 2 summarizes the floating ALC instructions 
with post-normalize option. 



Instruction 


Effect 


FMOV n, m 


FACn —FACm 


FPOSn.m 


1 FACn 1 — FACm 


FMNS n, m 


-|FACn| —FACm 


FNEG n, m 


- FACn —FACm 


FRND n, m 


rounded FACn — FACm 


FADD n, m 


FACn+FACm —FACm 


FSUBn.m 


FACm- FACn — FACm 



Table 2 - ALC Instructions with Post- 
Normalize Option 



FDIV n,m 
divides FACm by FACn and moves the result to FACm. 
The instruction 

FHLV n, m 

halves FACn and moves the result to FACm. 

Table 3 summarizes the ALC instructions that always 
post-normalize the result. 



Instruction 



FMPYn, m 
FDIVn, m 
FHLVn. m 



Effect 



FACm*FACn -FACm 
FACm/FACn —FACm 
FACn/2. — FACm 



Table 3 - ALC Instructions that Always 
Post- Normalize 

Floating ALC Instruction Examples 

Some examples of legal floating point ALC instructions 
are: 

FMPY 1,0 

FADD 0,1,FSGE 

FSUB# 1,0, FSLT 

FMOVU 3,0 

FMOV 0,0, FSLT 

FNEG 0,0,FSKP 



^T 



[_xj gives the maximum integer K such that K <c x , 
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Input/Output Instructions 

The use of I/O instructions requires the user to pro- 
vide two special routines. The first is an input routine 
which, when called, must return an ASCII input char- 
acter, right justified in AC0 with bit 8 = 0. The 
address of this routine must be stored by the user in 
location 040 of page Q. 

The second routine is an output routine which, when 
called, must accept an ASCII output character, right 
justified in AC0 with bit 8=0. The address of this 
routine must be stored in location 041 of page 0. All 
output messages to this routine will be terminated 
by a null (all zero) character. 

These I/O routines must be reentrant for the interpreter 
to be reentrant. (If the routines which interrupt and 
use the interpreter do not use I/O instructions, the 
user I/O routines need not be reentrant. ) 

The instruction 

FDFC n ;FLOATING POINT DECIMAL TO 
;FLOATlNG CONVERT 

will cause an ASCII character string in engineering 
notation to be converted to internal floating point form 
and loaded in FACn. The input characters must be 
provided by the user routine whose address is stored 
in location 40 8 of page 0. Numbers in the following 
form will be converted. 



.n[.]n... 



where : n is the decimal mantissa (the first seven non-zero 
digits will be converted and the remaining digits ig- 
nored) and each m is a digit of the decimal charact- 
eristic. 

The signs of the mantissa and characteristic are 
optional with the default assumed + . 

The break character is any character other than 

1. a decimal digit or 

2. an *'E", "+" , "-", or "." 

If the break character is a rubout (177), the entire 
string will be ignored and a new one must be given, 
i.e. , the conversion starts over. 

If the conversion results in a number less than 
16**-l*16**-63, the underflow flag will be set and 
true zero will replace FACn. If the conversion results 
in a number whose magnitude is greater than 





Input/Output Instructions (cont'd) 

(1-16**-6)*16**63, this latter magnitude will replace 
FACn and overflow will be set. As described on page 
2-1, "input conversion returns three additional words 
of information: conversion OK flag, the break char- 
acter, and decimal point seen flag. These may be 
examined and used as necessary. 

Examples of legal character strings are: 

1* 

1.* 
-1* 
+1* 

1E3* 

3.1415926* 

l.E+70* 

l.E-70* 

where * will be returned as the break character and 
conversion OK will be non-zero. 

Some illegal character strings are: 

A (break character will be A) 
+* (break character will be *) 
+. ! (break character will be !) 

Conversion OK will be zero in all these illegal cases. 

The instruction 

FDFCI n ; FDFC WITH JUDICATION 

will provide the user with an indication before the 
conversion begins. The ASCII character "F" followed 
by a null character will be passed to the output routine 
whose address is given in location 041. For example, 
if the user has provided for I/O from the teletype, 
the use of FDFCI will print "F" on the page copy every 
time an input is required. In all other respects it is 
identical to FDFC. 

The instruction 

FFDC n ; FLOATING POINT FLOATING TO 
;DECIMAL CONVERT 

will convert the number in FACn to an ASCII char- 
acter string in engineering notation. The output 
characters will be passed one at a time, right justi- 
fied in AC0, to a user routine whose address is stored 
in location 41 g of page 0. The output string will be of 
the form: 



/+>. nnnnnnnE/Amm 
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Input/Output Instructions (cont'd) 

where: each n represents a decimal digit of the man- 
tissa, 

each m represents a digit of the decimal 
characteristic. 

The string will be terminated by a null character. 

Special Instructions 

Two special instructions are defined which modify the 
index registers. 

The instruction 

FIC2 
causes AC2 to be incremented by two . 
The instruction 

FIC3 

causes AC3 to be incremented by two. These instruc- 
tions are useful for indexing through a table of float- 
ing point numbers. Use of F1C2 in indexing is shown in 
the second example in Chapter 5. 

A third special instruction provides a HALT feature 
within the interpretive mode. The instruction 

FHLT 

will cause the interpreter to HALT. Hardware AC" 1 
will contain the address of the FHLT instruction. The 
address lights will have no apparent relationship to the 
HALT, since the address is within the interpreter. The 
user may press CONTINUE to resume after this HALT. 

Illegal Instructions 

The proper encoding for all floating point instructions 
is given in Appendix B. The interpreter will HALT if 
an illegal instruction is encountered. Hardware AC0 
will contain the address where the illegal instruction 
was found. This HALT will occur if extended instruc- 
tions are used and only the Basic Interpreter is loaded, 
or on any bit configuration that cannot be decoded into a 
floating point instruction. The user cannot press 
CONTINUE to resume after this HALT. 



Chapter 2 - Basic Interpreter Use and Command Structure 

REQUIREMENTS FOR REENTRANCE 

A number of points regarding reentrance of the interpreter 
have been mentioned. This section explicitly defines 
the rules which must be followed by any routine which 
interrupts a base level routine and reenters the 
interpreter. 



1) 



2) 



3) 



All hardware accumulators, Carry, and 
page locations 006 and 007 must be 
saved, 

A new writable area address must be 
provided in location 007. 

An FINI must be issued after location 007 
has been set up (only necessary the first 
time). 



4) If I/O instructions are to be used, the 
user I/O routines must be reentrant. 
(Alternatively, locations 040 and 041 
must be saved and addresses provided to 
different I/O routines). 

5) Upon exit to the base level routine, the 
hardware accumulators, Carry and 
locations 006 and 007 must be restored. 
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FEATURES 



The Extended Interpreter pjw ides the instructions des- 
cribed in Chapter 2, plus a number of mathematical 
functions and "F" format output. If the Extended Inter- 

nvi ti^-r io nqoH lMrt //lQ^imrtl\ llr «^^ /,f .,.^JJ — Kid ~«-— -i 

^-.^^^^j. itJ tiS^vj, A W \u^^iinai/ vvuxuo Ul Wlllaulc CLutd^C 

must be provided by the user. 
MATHEMATICAL FUNCTIONS 



The mathematical functions are implemented using ALC 
instructions which always post-normalize. They permit 
the no load option as weii as the floating skip options. 
Appendix C provides a detailed description of the methods 



used to implement these functions as well as a discussion 
of their accuracy. The following is a general Ascription 
of each instruction. 



me instruction 

FALG n, m 

computes the natural logarithm of the contents of FACn 
and moves the result to FACm. If the argument is les"s 
than 0, the overflow flag is set and -(1-16**-6)*16**63 
is returned as the result. 

The instruction 

FATN n, m 

computes the arctangent of the contents of FACn and 
moves the result to FACm. The result is an angle 
expressed in radians in the range -tt/2<< arctan(x).< 

ir/2. 

The instruction 

FCOS n, m 

computes the cosine of the contents of FACn and moves 
the result to FACm. The argument is assumed to be an 
angle expressed in radians. If the argument is greater 
than 2**24, the overflow flag is set, and the result will 
be incorrect. 

The instruction 

FSIN n, m 

computes the sine of the contents of FACn and moves the 
result to FACm. The argument is assumed to be an 
angle expressed in radians. If the argument is greater 
than 2**24, the overflow flag is set,and the result will 
be incorrect. 



F LOATING POINT INTERPRETER 

Chapter 3 - Extended Floating Point Interpreter 

MATHEMATICAL FUNi TIONS (cont'd) 

The instruction 

FTAN n, m 

computes the tangent of the contents of FACn and moves 
the result to FACm. The argument is assumed to be an 
angle expressed in radians. If the argument is greater 
than 2**24, the overflow flag will be set,and the result 
will be incorrect. 

The instruction 

FEXP n, m 

computes e raised to the power contained in FACn 
and moves the result to FACm. If the argument is" 
less than -177.5, true is returned^and the under- 
flow flag is set. If the argument is greater than 
174.673, +(1-16**-6)*16**63 is returned.and the 
overflow flag is set. 

The instruction 

FSQR n, m 

computes the square root of the argument in FACn and 
moves the result to FACm. If the argument is less 
than 0, the underflow flag is set, and true is returned 
as the result. 

Table 4 summarizes the math functions. 



Instruction 


Effect 


FALG n, m 


ln(FACn) — FACm 


FATN n, m 


arctan (FACn) —FACm 


FCOS n, m 


cosine (FACn) —FACm 


FSIN n, m 


sine (FACn) —FACm 


FTAN n, m 


tangent (FACn) —FACm 


FEXP n, m 


e**(FACn) -FACm 


FSQR n, m 


(FACn)**l/2 -FACm 



Table 4 - Math Functions 
'F" Format Conversion 



The Basic Interpreter provides floating point to decimal 
conversion with "E" format output. The extended ver- 
sion provides "F" format output as well. 

The instruction 

FFDCF n ;FFDC WITH "F" FORMAT OUTPUT 
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Chapter 3 - Extended Floating Point Interpreter 

"F" Format Conversion (Cont'd) 

will convert the floating point number in FACn to decimal 
and output a character string in "F" format via the user 
routine whose address is stored in location 041. The 
output will be of the form: 

AAA&-1 nnn , nnnn 
C \*- D ■ - 



W 



The width of the field (including sign and decimal point), 
W, and the number of places to be given after the dec- 
imal point, D, must be set up in the writable area before 
FFDCF is given. The displacements of these words are 

W = 121 8 
D = 122 g 

They can be accessed in. a manner similar to the flags 
described on page 2-1. Two conditions will cause the 
overflow flag to be set and no conversion to be performed. 



"F" Format Conversion (cont'd) 

1. W > 32 (entire width of field limited to 
32 characters. ) 

2. W < D + 2 (W must be 2 greater than D 
to provide for sign and decimal point, ) 

If W is not large enough to accommodate the number, 
significant digits will be lost. 

Assume W has been set to 12 (decimal) and D to 6. 
Examples of "F" format versus "E" format outputs are: 



'E" Format 



+. 1374600E+02 
-.7968433E-03 
+. 1000000E+04 
-.1 000000 E+05 
+. 1000000E+06 
-.3500000E-06 
+.4713279E-01 



"F" Format 



+13.746000 

-.000796 
+1000.000000 

10000.000000 (note sign lost) 
00000. 000000 (all lost) 

-.000000 (significance lost) 

+.047132 
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The Extended Relocatable Floating Point Interpreter is 
identical to the Extended Floating Point Interpreter, ex- 
cept for the following: 

1 ) No absolute locations are used. 

2) The interpreter requires 3 zero relocatable lo- 
cations* and approximately 3500g normal reloca- 
table locations. 

3) Within any program that calls or initializes the 
interpreter, the appropriate normal external 
mnemonics must be declared. These are* 

FENT -FLOATING INTERPRETER ENTER 
FINT -FLOATING INTERPRETER 

•INITIALIZE 

FENT replaces the FRTP_ command , used to call 
the absolute interpreter. 

FINT replaces the FINI command, used to 
initialize the absolute interpreter. 

For example: 

.EXTN FENT, FINT 

.NREL 

FINT ; INITIALIZE 

FENT -ENTER INTERPRETER 

FLDA 0, FONE ; FLOATING INSTRUCTIONS 



Chapter 4 - Extended Relocatable Floating Point Interpreter 

5) If I/O instructions are to be used, two more 
ZREL locations must be defined and declared 
as entries with the labels: 

GETC 
PUTC 

Location GETC must contain a pointer to the 
input character routine. GETC replaces 
location 040 of the absolute interpreter. 

Location PUTC must contain a pointer to the 
output character routine. PUTC replaces 
location 041 of the absolute interpreter. 



4) One ZREL location must be defined and declared 
as an entry (. ENT) with the label 

WSA 

Location WSA must contain a pointer to the 
writable area of the interpreter. WSA replaces 
location 0p7 of the absolute interpreter. (The 
writable area must be 1W\Q contiguous locations 
as in the absolute interpreter. ) 



6) Requirements for Reentrance 

a) All hardware accumulators, Carry, and 
tuc contents of WSA and ITMF must be 
saved. 

b^ A new writable ar&a arldrooo t™,^«- k~ 

provided in WSA. 

c) A FINI must be issued after location 
WSA has been set up. (This is necessary 
only once per writable area. ) 

d) If I/O instructions are to be used, the 
user I/O routines must be reentrant (or 
locations GETC and PUTC must be 
saved and new addresses provided to 
other I/O routines. ) 

e) Upon exit, all of the information saved 
above must be restored. 

The Extended Relocatable Floating Point Interpreter can 
be loaded for stand-alone use or for use under the Disk 
Operating System. The relocatable version of the 
interpreter is sent in relocatable binary format, Tape 
#089-000046. 



Three ZREL locations are required in the interpreter, and three are provided by the user program (WSA, 
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Chapter 5 - Sample Programs 



SQUARE ROOT NEWTON ITERATION 



The following routine is an example of a floating point subroutine that perform 
a square root Newton iteration on the trial guess in FAC0 given the argument 
FAC1. 



THIS ROUTINE PERFORMS A SQUARE ROOT NEWTON ITERATION 
FAC1 CONTAINS THE ARGUMENT AND IS NOT DESTROYED 
FAC0 CONTAINS THE PREVIOUS GUESSTIMATE 
COMPUTES (FAC0 + (FAC1/FAC0)) /2. 
CALLING SEQUENCE 

FTSR NSR 

RETURN : RESULT IN FAC0 



NSR: 



FMOV 1,2 
FDIV 0, 2 
FADD 2,0 
FHLV 0, 
FJMP0.3 



SAVE ARGUMENT 
FAC1/FAC0 

(FAC0+FACl/FAC0)/2. 
RETURN 



POLYNOMIAL EVALUATION 

This polynomial evaluation routine requires AC2 to point to the first word of a 
table of floating point coefficients, ordered high order coefficient down and 
terminated by true 0. For example: 



AC2 



Cn 






Cn-1 






Cn-2 






* 




/ 



Word 



1 
2 
3 
4 





/ 




■ 




CI 




C0 








a 




- 


Y> 




- 


- - - - - 


- - 



2n-f# 
2n+l 
2n+2 

2n+3 



The routine uses Horner's method for evaluation, i.e. 

f(x) = (. . . ( ( (x + C n ) x + Cn-i) x + Cn-2) x . . .+ Cj) x + C# 
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Chapter 5 - Sample Programs 

POLYNOMIAL EVALUATION (Continued) 



POLYNOMIAL EVALUATION 

FAC2 CONTAINS ARGUMENT (X) 

AC2 POINTS TO COEFFICIENT LIST TERMINATED BY 0, 

AND ORDERED HIGH TO LOW. 
RESULT RETURNED IN FAC0 
FAC0, FAC1 DESTROYED 
CALLING SEQUENCE 

FJSR FPLY 

RETURN 



FZRO: 














FLPY: 


FLDA 


0, FZRO 


FLPY1 : 


FLDA 


1,0,2 




FMOV 


l.l.FSNR 




FJMP 


0,3 




FMPY 


2,0 




FADD 


1.0 




FIC2 






FJMP 


FPLY1 



;CLEAR RESULT 
;GET COEFFICIENT 

RETURN IF ZERO 

SUM * ARGUMENT 

SUM * ARG. + COEF. 

BUMP POINTER TO NEXT COEF. 
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Appendix A - Floating Point Instructions 



Standard Instructions 



# 

FSGE 

FSGT 

FSKP 

FSLE 

FSLT 

FSNR 

FSZR 

U 



FADD 


Floating Add 


FDFC 


Floating Decimal to Floating Convert 


FDFCI 


FDFC with Indication ~ 


FDIV 


Floating Divide 


FDSZ 


Floating Decrement and Skip if Zero 


FETR 


Floating Mode Enter 


FEXT 


Floating Mode Exit 


FFDC 


Floating Floating to Decimal Convert 


FFIX 


Floating to Fixed 


FFLO 


Fixed to Floating 


FHLT 


Floating Halt 


FHLV 


Floating Halve 


FIC2 


Floating Increment AC2 


FIC3 


Floating Increment AC3 


FINI 


Floating Initialize 


FISZ 


Floating Jncrement and Skip if Zero 


Fjivip 


Floating Jump ~~ 


FJSR 


Floating Jump to Subroutine 


FLD3 


Floating Load AC3 


FLDA 


Floating Load Floating Accumulator 


FMOV 


Floating Move 


FMNS 


Floating Move Minus 


FMPY 


Floating Multi£l^ 


FNEG 


Floating Negate 


FPOS 


Floating Move Positive 


FRND 


Floating Round 


FST3 


FloatingjStore AC3 


FSTA 


Floating _Store Floating Accumulator 


FSUB 


Floating Subtract 


Extended Instructions 




FALG 


Floating Natural Logarithm 


FATN 


Floating Arctangent 


FCOS 


Floating Cosine 


FEXP 


Floating Exponential 


FFDCF 


Floating Floating to Decimal Convert \ 


FSIN 


Floating Sine 


FSQR 


Floating Square Root 


FTAN 


Floating Tangent 


Floating Point Options 





No Load 

Floating Skip on Greater Than or Equal 

Floating Skip on Greater Than 

Floating Skip 

Floating Skip on Less Than or Equal 

Floating Skip on Less Than 

Floating Skip on Non Zero Result 

Floating Skip on Zero Result 

Unnormalize (no post-normalization) 
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Appendix B - Floating Point Instruction Encoding 



MEMORY REFERENCE WITHOUT ACCUMULATOR 

INDIRECT 






FUNCTION 


V 


INDEX 


DISPLACEMENT 





1 2 1 1 4 


s 


. 1 , 


• V f 10 II 12 1 13 14 IS 



0000 FJMP 1100 FFL0 

0001 FJSR 1101 FLD3 

0010 FISZ 1110 FST3 

0011 FDSZ 1111 FFIX 

MEMORY REFERENCE WITH ACCUMULATOR 

INDIRECT 



FUNC- 
TION 



FAC 
ADDRESS 



INDEX 



DISPLACEMENT 



I 2 3 I 4 S 4 I 7 ■ 9 I 10 II 12 | 13 |4 IS 



01 FLDA 
10FSTA 



SPECIAL 



1 



FUNC- 
TION 



00000000000 



• I I 1 J I 4 S 4 17 t 4 I 10 II 12 I 13 14 15 



00FEXT 
01 FIC2 
01 FIC3 
11FHLT 









CONVERSION 




1 


FUNC- 
TION 


FAC 
ADDRESS 


0000000000 


TYPE 





1 2 


3 .4 


S * 1 7 1 9 1 It II 12 1 13 14 


15 



01 FDFC 
10 FFDC 



FDFC 
II 

FFDC 
1 F 



ARITHMETIC (OPTIONAL NORMALIZATION) 



NO 

NORMALIZE 



FAC 
SOURCE 
ADDRESS 



FAC 

DEST. 

ADDRESS 



FUNCTION 







I 2 3 14 5 4 17 > 1 I 10 II 12 



SKIP 



001 FNEG 

010 FM0V 

011 FP0S 

100 FMNS 

101 FSUB 

110 FADD 

111 FRND 



001 FSGT 

010 FSLT 

011 FSNR 

100 FSZR 

101 FSGE 

110 FSLE 

111 FSKP 



ARITHMETIC (ALWAYS NORMALIZED) 



1 


FAC 
SOURCE 
ADDRESS 


FAC 

DEST. 

ADDRESS 




' 1 ' 


FUNCTION 


NO 
LD. 


SKIP 



.10 II 12 



13 M 1$ 



0001 FALG 


001 


FSGT 


0010 FATN 


010 


FSLT 


0011 FCOS 


011 


FSNR 


0100 FMPY 


100 


FSZR 


0101 FSIN 


101 


FSGE 


0110 FTAN 


no 


FSLE 


1000 FDIV 


in 


FSKP 


1001 FEXP 






1010 FSQR 






1100 FHLV 
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This appendix describes in some detail the function, 
method, and accuracy of the mathematical routines 
supplied with the Extended Interpreter, The accuracy 
of the routines is influenced by two factors: 

1. The accuracy of the argument 

2. The accuracy of the algorithm 

These two factors will be mentioned for each routine. 
The accuracy of the algorithm itself assumes an 
argument that is exact, i. e. , no argument error. 

The relative and absolute errors of a function routine 
are defined as follows: 

Let f(x) = the true value of the function at x 

^ct gyX/ = tuC resuit returned by uie function 
routine, given x 

Now the absolute error- <->f th& rponit ;o 



ABS (f(x) - g(x)) 
and the relative error of the result is 

ABS ((f(x) - g(x))/f(x)) 
FLOATING POINT ARCTANGENT 
Function: 



To calculate the arctangent of x, where x is a floating 
point number, and return an angle in radians in the 
range -tt/2 <c arctan (x) < ir/2 . 

Method : 

The range of x is immediately reduced to 

^ x < 1 
by means of the identities 

arctan (-abs(x)) = - arctan(abs(x)) 
and if abs(x) > 1, 

arctan (abs(x)) = T/2-arctan(l/abs(x)) 
For x > tan ( tt/12) the range is reduced to 

tan (-T/12)^y^tan(7r/12) 



Appendix C - Extended Interpreter Mathematical Routines 
FLOATING POINT ARCTANGENT (Cont d) 

u y jiiGuJIS Oj. Li.iC IGCIitity 

arctan(x) = arctan(x0) + arctan [(x-x0)/(M-x*x0)]. 
For x0 = 1/3**0. 5, we obtain 

arctan (x) = tt /6 + arctan [((x*3**0. 5)-l)/(x+3**0. 5)] 
= 7r/6 + arctan (y) 

where y satisfies the range given above. The arctan 
is computed using the first four terms of a polynomial 
approximation of the form: 

n 
arctan (x) = x * T~ C *x**(2*i) 
i=0 



Accuracy : 



Argument nrror 



If x is the argument, the absolute error of the result 
is approximately 

£/(l+x**2) 

where £ is the absolute error in x. Thus for small 
values of x, the errors are almost equal, while as 
x becomes larger, the effect of the argument error 
decreases. 

Maximum Relative Error 

For the range 

-tan (V12) ^c x x < tan(7/12), 

the maximum relative error is approximately 
10**- 7. 6 . 

FLOATING POINT EXPONENTIAL 
Function : 

To calculate e to the power x, where x is a floating 
point number. 

Method : 

If x < -177. 5, return as the result and set underflow 
flag. 

If x > 174.673, return the largest positive number as 
the result and set overflow flag. 
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Appendix C - Extended Interpreter Mathematical Routines 

FLOATING POINT EXPONENTIAL (cont d) 

Otherwise, we use the equality 

e ** x = 2**(x*log2(e)) 

Let x*log2(e)= m*16**c 

Further, let q = L.m*16**cJ and f = (m*16**c)-q. 
Therefore ^ f < 1.0 . 

Now2**(f +q)= 2**f*2**q 

Let us compute 2 * *f . 

The range of f can be further reduced if we let 

g = f - 0. 5 where -0. 5 < g < 0. 5 . 

We compute 2**g directly if g >, , otherwise we compute 
2**g as 

l/(2**abs(g)) 

2 **g (where ,£g <0. 5) is computed using the first five 
terms of a polynomial approximation of the form: 

n 
2**g = £ C.*g**i 

i=0 

Now f = g + 0. 5 

Therefore, 2**f = 2**g*2**0.5 

But the answer is 2**f*2**q 

q is an integer and we let q = 4i + j 

Now 2 **f*2**q = 2**f*16**i*2**j 

The characteristic of f is added to i to obtain 2**f*16**i. 
This result is shifted j positions left if i > or j 
positions right if i <0. The result is, of course, e**x. 



Accuracy: 



Argument Error 



The relative error of the result is approximately equal 
to the absolute error of the argument. Thus for large 
values of x, substantial relative errors in the results 
can occur. 

Maximum Relative Error 

For the range < x < 0. 5 the maximum relative error of 
2 **x is approximately 10**- 7. 0. 



FLOATING POINT NATURAL LOGARITHM 

Function : 

To calculate the natural logarithm of x, where x is 
a floating point number. 

Method : 

If x < 0, overflow flag is set and minus the largest 
floating point number is returned as the value. 
Otherwise, let x = m*l6**c. By means of a binary 
normalization, the range of m' is reduced to 

l/2<cm'<l , 

and x = m' * 16 **p*2**(-q) where q = number of left 
shifts required to normalize (0 »< q £ 3 ). 

Now for 1/2 £ m' < 1/2 * *0. 5 

let a = 1/2, b = 1 

for 1/2**0. 5 <:m'<l 

leta = 1, b= 

Define y = (m' - a)/(m'+ a) 

then m' = a* (y+l)/(-y+l) 

Nowx= 16**P*2**(-q)*a*(l+y)/(l-y) 
= 2**(4p-q-b)*(l+y)/(l-y) 

Using ln(x) = ln(2) * log2(x) we obtain 

ln(x)= ln(2)* [(4p-q-b>Hog2((l+y)/(l-y))] 

ln(x) = ln(2) * (4p-q-b) +ln((l+y)/(l-y )) 

From the above, we can determine that 

2**-0.5 < (l +y )/(l-y) $ 2**0.5 

The In ((l+y)/(l-y)) is determined for the above range using 
the first three terms of a polynomial approximation of the 
form: 

n 
ln( z) =? z * Y, <Vz**(2*i) 

i'=0 

where z = (l+y)/(l-y). 
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FLOATING POINT NATURAL LOGARITHM (con t) 

Accuracy : 

Argument Error 

The absolute error in the result is approximately equal 
to the relative error in the argument. Therefore, an 
argument close to 1 can give a large error since the 
function at this value is quite small. 

Maximum Relative Error 

For the range 

l/2**-0. 5 ,<(l+x)/(l-x) «2**0. 5 

the maximum absolute error of ln(x) is approximately 
in**-7 6 



FLOATING POINT SINE AND COSINE 
Function : 

To calculate sin(x) or cos(x), where x is the floating 
point angle in radians. 

Method : 

Compute p = abs(x)*4/ir 

Let q = j_pj , f = p-q where ^ f < 1 

Now q represents the half -quadrant in which the abs(x) 
falls. Using the following equalities. 

sin(x) = -sin(-x) 
cos(x) = cos(-x) 
cos(x) = sin(x+V2) 

we define 

ql = q if sin is required and x ^ 

ql = q + 2 if cos is required 

ql = q + 4 if sin is required and x < 

Then for all values of x, the computation has been 
reduced to 

sin(7y'4*(q 1 +f))= sin(t) 

Since sine (and cosine) are periodic in 2*74 we take q-. = 
q^ mod 8. Using the further equality 

sine(7j/4+x) = cosine(7^4-x) 



Appendix C - Extended Interpreter Mathematical Routines 
FLOATING POINT SINE AND COSINE . (cont d) 
we finally can produce the table 



qi 


sin(t) 





sin(V4*f) 


1 


cos(V4*(l-f)) 


2 


cos(7/4*f) 


3 


sin(7^4* (1-f)) 


4 


-sin(V4*f) 


D 


-cos(V4*(l-f)) 





-COS(1f4 t) 


7 


-sin(V4*(l-f)) 



In all cases, the argument has been reduced to the 
range jet < it /4 . The sine is computed using the 
first four terms of a polynomial approximation of 



uic 10x111: 






i = 



r 1 *„**/1 *:\ 



The cosine is computed using the first four terms of 
a polynomial approximation of the form: 

n 
cos(x) = Y. C i *x**(2*i) 

i = 



Accuracy : 



Argument Error 



The absolute error of the result is approximately equal 
to the absolute error in the argument. Thus, the 
larger the argument, the larger the absolute error 
of the result. 

Maximum Relative Error 

For the range .< x < V4 the maximum relative 
error for sin(x) and cos(x) is approximately 10**-7. 4 . 

FLOATING POINT SQUARE ROOT 

Function : 

To calculate the square root of a floating point 
number, x. 

Method : 

Let x= m*16**c 



C-3 



FLOATING POINT INTERPRETER 



Appendix C - Extended Interpreter Mathematical Routines 

FLOATING POINT SQUARE ROOT (cont'd) 

If m < 0, set underflow flag and return 0. as the result. 
If m = 0, return 0. as the result. 

Otherwise, let c = 2*p + q where p is an integer and 
q = or 1. Now if q = 0, we have x = m*16**2p 
and x**l/2 = m**l/2*16**p. 

If q = 1, we have x = m*16**(2p+l) 

= m*16**(2p+2)/16 

and x**l/2 = (m**l/2)/4*16**(p+l) 

Therefore, the characteristic of the result is p + q, 
and the problem has been reduced to finding a suitable 
first guess tor m**l/2 if q = or (m**l/2)/4 if q = 1. 

An initial guess is taken in the hyperbolic form 
y0 = a+b/(c+m) 



where for q = 0, 



and for q = 1, 



a = 1.80713 
b= -1.57727 
c = 0.954182 

a= 0.428795 
b = -0. 3430368 
c = 0. 877552 



The initial guess is now 

yl = y0*16**(p+q) 

Two Newton iterations give us the result. 

y2 = (yl+x/yl)/2 
x**l/2 = y3 = (y2+x/y2)/2 



Accuracy: 



Argument Error 



The relative error of the result is approximately half 
the relative error in the argument. 

Maximum Relative Error 



The maximum relative error for x**0. 5 is approximately 
10(**-6.0. 



•LOATING POINT TANGENT 



Function: 



To calculate the tangent of x, where x is the floating 
point angle in radians. 

Method : 

Compute p = abs(x)*4/Tr 

Let q = i p I , f=p-q where ,< f < 1 

Take qj = q mod 4 and in a manner similar to sine- 
cosine we can obtain the table below: 



_2L 



tan(x) 



tan (ir/4*f) 

cot (7T/4*(l-f)) 

-cot (7r/4*f) 
-tan (tt/4*(1-i)) 



In all cases, the argument has been reduced to the 
range <; arg < 7r/4 . 

The tangent is computed using the first six terms of 
a polynomial approximation of the form: 

n 
tan(x) = x*£ C^x**^*!) 

i = 



Accuracy : 



Argument Error 



The absolute error of the result is approximately 
equal to 

£* (l+tan(x)**2) 

where £ is the absolute error of the argument. Thus 
if x is near an odd multiple of w/2, an argument error 
will produce a larger absolute error in the result. 

Maximum Relative Error 



For the range 

4: x < tt/4 

the maximum relative error of tan(x) is approximately 
10**-6.6 . 



C-4 



CHANGES FROM REVISION 3 TO REVISION 4 OF THE FLOATING POINT 
INTERPRETER USER'S MANUAL 



Page 
1-1 



Lower bound of the decimal range of a floating point number changed 
to2.4xl0' 78 . 



Spelling of during corrected. 

2-1 Second note at bottom of page corrected. 

2-4 The action taken by FMOV 0, 1 is corrected to read ". . . FMOV 0, 1 

moves FACO to FAC1 ..." 

2-7 Sentence on indexing example corrected. 
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ARCTANGENT C-l, 3-1 

ARITHMETIC and logical instructions 2-4 



BASIC Interpreter Chapter 2 
BREAK character flag 2-1 



CHARACTER string conversion 2-6 
CONVERSION 

character string 2-6 

E format 3-2 

F format 3-1 

OK flag 2-1 
COSINE C-3, 3-1 



DECIMAL point seen flag 2-1 



E Format conversion 3-2 

EFFECTIVE address 2-3 

ENCODING of instructions App. B 

EXPONENTIAL C-l, 3-1 

EXTENDED Interpreter (binary) Chapter 3 

EXTENDED Interpreter (relocatable) Chapter 4 



F Format conversion 3-1 
FENT mnemonic 4-1 
FINT mnemonic 4-1 
FLAGS 
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overflow/ underflow 2-1 
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FSGT skip option 2-4 
FSKP skip option 2-4 
FSLE skip option 2-4 
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FSNR skip option 2-4 
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INITIALIZATION 2-2 
INPUT/OUTPUT Instructions 2-6 
INSTRUCTIONS 
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illegal 2-7 
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input/output 2-6 
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FADD ' 2-5 
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FCOS 3-1 
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FDIV 2-5 



FDSZ 
FETR 



2-3 

2-2, 2-3 
3-1 



GETC (get character) 4-1 



HALT 2-7 



INDEX register 

modification 2-7 

use in instructions 2-1 



FEXT 2-2 
FFDC 2-5 
FFDCF3-1 
FFIX 2-3 
FFLO 2-3 
FHLT 2-7 
FHLV 2-5 
FIC2 2-7 
FIC3 2-7 
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FISZ 2-3 
FJMP 2-3 
FJSR 2-3 
FLD3 2-3 
FLDA 2-1 
FMOV 2-4 
FMNS 2-4 
FMPY 2-5 
FNEG 2-5 
FPOS 2-4 
FRND 2-5 
FSIN 3-1 
FSQR 3-1 
FST3 2-3 
FSTA 2-3 
FSUB 2-5 
FTAN 3-1 
INTERPRETER 

Basic 2-1, Chapter 2 
Definition of i 
entering the 2-2 
exiting the 2-2 
Extended 2-1, Chapter 3 
Extended relocatable Chapter 4 
general structure 2-1 
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linkage enabling 1-1 
requirements for 2-7, 4-1 
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SAMPLE programs Chapter 5 

SINE C-3, 3-1 

SKIP options 2-4, 3-1 
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initialization 1-1 

interpreter 1-1 
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